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As the theme for the future L3 Cosmic Vision mission, ESA has recently chosen the ‘Gravita¬ 
tional Wave Universe’. Within this call, a mission concept called eLISA has been proposed. This 
observatory has a current initial configuration consisting of 4 laser links between the three satellites, 
which are separated by a distance of one million kilometers, constructing a single channel Michelson 
interferometer. However, the final configuration for the observatory will not be fixed until the end 
of this decade. With this in mind, we investigate the effect of different eLISA-like configurations 
on massive black hole detections. This work compares the results of a Bayesian inference study of 
120 massive black hole binaries out to a redshift of « ~ 13 for a 10®m arm-length eLISA with four 
and six links, as well as a 2 x 10®m arm-length observatory with four links. We demonstrate that 
the original eLISA configuration should allow us to recover the luminosity distance of the source 
with an error of less than 10% out to a redshift of z ~ 4, and a sky error box of Afl < 10^ deg^ out 
to 2 ~ 0.1. In contrast, both alternative configurations suggest that we should be able to conduct 
the same parameter recovery with errors of less than 10% in luminosity distance out to 2 12 

and AU < 10^ deg^ out to 2 ~ 0.4. Using the information from these studies, we also infer that 
if we were able to construct a 2Gm, 6-link detector, the above values would shift to 2 20 for 

luminosity distance and 2 ~ 0.9 for sky error. While the final configuration will also be dependent 
on both technological and financial considerations, our study suggests that increasing the size of a 
two arm detector is a viable alternative to the inclusion of a third arm in a smaller detector. More 
importantly, this work further suggests no clear scientific loss between either choice. 


PACS numbers: 04.30.Tv,95.85.Sz,98.90.Es,98.62.Js 


I. INTRODUCTION 

The European Space Agency (ESA) has recently cho¬ 
sen the theme of the “Gravitational Wave Universe” for 
the Cosmic Vision L3 mission selection. Within this 
program, a gravitational wave (GW) observatory called 
eLISA has been proposed [1, 2]. This observatory is com¬ 
posed of three space-craft forming an equilateral triangle. 
The proposed constellation is a single-channel Michelson 
interferometer, made up of one mother, and two daugh¬ 
ter spacecraft separated by 10® kms. It is intended that 
this observatory will function in the frequency band of 
10“® < f /Hz < 1. Such a mission should be capable 
of observing GWs from stellar mass compact binaries in 
our own galaxy, the inspiral of stellar mass objects into 
supermassive black holes (the so-called extreme mass ra¬ 
tio inspirals or EMRIs) out to a redshift of 2 ^ 0.7, 
and possibly even cosmological defects such as cosmic 
(super)strings or a stochastic cosmological background. 
However, out of all the possible sources, the most abun¬ 
dant, and brightest, should be the merger of supermas¬ 
sive black hole binaries (SMBHBs) out to a large redshift. 

While there is evidence to suggest the existence of su¬ 
permassive black holes at the center of each galaxy, and 
also evidence to suggest an intricate connection between 
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the mass and evolution of the SMBH and the evolution 
of the host galaxy [3-5], our knowledge of the formation 
and evolution of massive black holes is slim. It is possible 
that the first black holes, either formed as remnants from 
the first short-lived, low metallicity, stars (called Pop III 
stars) at redshifts of 2 ^ 20 [6-13], or from the direct 
collapse of protogalactic disks at redshifts of 2 < 12 [14- 
17]. In the first scenario, black holes with masses of a 
few tens to a few thousand solar masses would have been 
produced, while in the second scenario, the initial black 
holes would have had masses of 10® — 10® Mq. 

It is expected that with eLISA, it should be possible 
to measure the system parameters with sufficient accu¬ 
racy that one could investigate the models of BH seed 
formation and evolution [18, 19], as well as possible de¬ 
viations from General Relativity [20-22]. However, while 
the theme for the L3 mission has been fixed, the final 
composition of the mission configuration will not be fixed 
until the end of this decade. This leaves time for the 
study of alternative mission configurations, and an inves¬ 
tigation of the impact it has on the science of SMBHBs. 

In a recent work, it was demonstrated that to unlock 
the full science potential of a space-based GW observa¬ 
tory, it is necessary to carry out a full Bayesian analysis 
when estimating the parameters of the binary [23] (here¬ 
after referred to as PC15). This study demonstrated that 
it is possible, with the proposed eLISA configuration of 
detecting SMBHBs out to a redshift of 2 ~ 13, just using 
a post-Newtonian (PN) inspiral waveform. The goal of 
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this work is to make a direct comparison between these 
recent results for eLISA, and some possible alternative 
mission configurations. 

A. Outline of the paper 

The paper is structured as follows. In Sec II we out¬ 
line the composition of an eLISA-like observatory, and 
define the instrumental noise model to be used in this 
study. In Sec III we define the post-Newtonian higher 
harmonic corrected inspiral waveform and the detector 
response to an impinging GW. We then use this wave¬ 
form and an astrophysically motivated binary black hole 
population to calculate the detection horizon for each of 
the configurations. In Sec IV we present a comparison 
between the original eLISA configuration and the two al¬ 
ternatives. Using a regression analysis, we also infer the 
scientific performance of a larger three arm detector. 



FIG. 1: A comparison of instrumental noise model for a IGm 
armlength, 4 and 6-link configuration (G1,C2) and a 2Gm 
armlength, 4-link, eLISA configuration (C3), assuming a con¬ 
stant laser power and telescope size. 


II. THE eLISA-LIKE DETECTOR RESPONSE 

The “Gravitational Wave Universe” has recently been 
chosen by the European Space Agency (ESA) as the 
theme for the L3 mission within the Gosmic Vision pro¬ 
gram. During the initial phase of the L3 competition, 
a concept for a GW observatory called eLISA/NGO [? 
] was proposed by the scientific community. This re¬ 
structured version of the LISA mission concept was again 
composed of three spacecraft in an equilateral triangular 
configuration, inclined at an angle of 60° to the plane of 
the ecliptic, with each spacecraft following heliocentric 
ballistic orbits. This causes the constellation to natu¬ 
rally cartwheel once per orbit. This motion induces a 
Doppler phase shift into the GW signal and is the princi¬ 
pal component in the sky resolution of a GW source. For 
the L3 mission, the proposed observatory was composed 
of a single mother and two daughter spacecraft, each sep¬ 
arated by 10® km. This configuration corresponds to a 
4-link, single channel, Michelson interferometer. 

As stated in the Introduction, the final mission config¬ 
uration for the L3 mission will not be fixed until around 
2020. As SMBHBs are the primary driving source for the 
mission design, in this work we will investigate the effects 
of two variants of the eLISA mission on the detection and 
parameter estimation for these sources. The first config¬ 
uration is identical in size to the eLISA proposal, except 
this time instead of the one mother-two daughter configu¬ 
ration, we have three mother spacecraft. This introduces 
another two laser links and corresponds to a two-channel 
Michelson interferometer with IGm armlengths. The sec¬ 
ond configuration is again a four-link, single channel in¬ 
terferometer, but this time with the armlengths doubled 
to 2Gm. In the rest of this article, for brevity, we will 
refer to the IGm, 4-link, eLISA as configuration Gl, the 
IGm, 6-link, configuration as G2, and the 2Gm, 4-link, 
observatory as configuration G3. 


Assuming that we keep the telescope size {D = 20 cm) 
and laser power [P = 2 W) constant, we can analytically 
describe the noise spectral density of the detector using 
the expression 
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where L is the arm-length of the particular eLISA con¬ 
figuration, S'P°®(/) = 1.21 X 10“^^m^/Hz is the position 
noise, S'i)°°(/) = 9 x 10“®®m^/(s'^Hz) is the acceleration 
noise, = 6.28 x 10“^®m^/Hz is a frequency indepen¬ 
dent fixed level noise in the detector and /* = l/(27rL) 
is the mean transfer frequency of the constellation. To 
count for possible problems in achieving design sensitivity 
at low frequencies, we have introduce a red noise compo¬ 
nent at frequencies less than 10“"^ Hz. 

In Fig (1) we plot the power spectral density of the 
instrumental noise as a function of frequency for the 1 
and 2Gm configurations. For the 2Gm detector, at lower 
frequencies we gain a factor of 2 in sensitivity as we can 
see from Eq (1) that \JSn{f) oc 1/L in power spectral 
density. As expected, by making the entire constella¬ 
tion larger, the bucket of the noise curve also moves to 
lower frequencies, making the 2Gm configuration sensi¬ 
tive to more massive black hole binaries. At high fre¬ 
quency there is little observable difference between the 
two noise curves as we are dominated at these frequen¬ 
cies by the photon shot noise from the laser. 
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III. THE BLACK HOLE DETECTION HORIZON 


A useful quantity for measuring the difference between 
observatory configurations is the detection horizon limit 
for SMBHBs (i.e. the maximum redshift, as a function of 
total redshifted mass, that a source can be detected at, 
given a detection threshold signal-to-noise ratio (SNR)). 
In PC 15 it was demonstrated that the Cl configuration 
had a maximum detection horizon oi z ^ 13 for systems 
with redshifted total masses of m{z) ^ 10^ Mq, given a 
SNR detection threshold of p = 10 [22], where SNR for 
a single channel interferometer is defined by 


^ {s\h) 


( 2 ) 


The angular brackets correspond to the usual noise 
weighted inner products 
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where s(/) and h{f) are the Fourier transforms of the 
time domain waveforms s(t) and h{t). 

In this work, we assume that the low frequency approx¬ 
imation (LFA) [24] is viable, and the GW wavelength is 
greater than the arm length of the detector. In this case, 
the strain of the detector to an impinging GW is 

h,it) = h+{OF+it) + h^iOF:^{t), (4) 

where i = { 1 , 2 } defines the channel number, 
the two polarizations of the CW, F^"^(t) are the beam 
pattern functions of the detector (which will be defined 
later) and f (t) is a phase shifted time parameter 

^(t) =t — sin 0 cos {a{t) — (j)), (5) 


where t is the time in the solar system barycenter, = 
\A\Jjc is the radial distance to the detector guiding cen¬ 
ter, c is the speed of light, ( 0 , cj)) are the position angles 
of the source in the sky, a(t) = 2'Kfmt + k, fm = 1/year 
is the constellation modulation frequency and k gives the 
initial ecliptic longitude of the guiding center. 

The beam pattern functions of the detector F+’^(t) 
are given in the low frequency approximation by 

F+{t\i/,e,(l)) = i [cos(2^/;)i:)+(t;'i/’,6',(/>, Ai) 

- sin(2'0)T>'^(t;-!/',6',<(), Ai)] , (6) 

+ cos{2/;)D^{t;ip,0,(l),X,)] , (7) 


where '0 is the polarization angle of the wave and we 
take Ai = 0 for the single channel configuration, or 
Xi = ( 0 , 7 r/ 4 ) for the two-channel configuration. Explicit 



FIG. 2: Massive black hole redshift horizon distance (exclud¬ 
ing merger and ringdown) as a function of redshifted total 
mass for configurations Cl (solid blue line), C2 (dashed red 
line) and C3 (dot-dashed orange line). Each curve assumes 
an inspiral-only detection threshold of pth = (10, 9, 8) respec¬ 
tively. 


expression for the detector pattern functions can 

be found in [25]. 

As with PC15, we will use the post-Newtonian (PN) 
inspiral waveform to 2PN order in phase and frequency, 
including higher harmonic corrections (HHCs) up to 2- 
PN order [26] 


l-l-.x — 
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where m = mi -I- m 2 is the total mass of the binary, 77 = 
mim 2 /m^ is the reduced mass ratio, x = (Gmuj/ 
is the invariant PN velocity, where to = d^orb/dt is the 
2 PN order circular orbital frequency as a function of 

(n) 

the orbital phase, the functions H/_ contain the PN 
corrections to the amplitude and the phase harmonics 
and Dl is the luminosity distance of the source. The 
luminosity distance, Dl{z), is defined within a ACDM 
model by 

D _ n 4- 

Ho Jo A/nfl:(l -I- z')‘^ + nM(l + z')^ + Ha ’ 

( 9 ) 

where the concurrent PLANCK values of ll/j = 4.9 x 
10“^, = 0.3086 and = 0.6914 and a Hubble’s 

constant of Hq = 67 km/s/Mpc are used [27]. 

The first step of our analysis is now to set a detection 
threshold for configurations C2 and C3. To do this, we 
carry out null tests where the output of our detector is as¬ 
sumed to be pure Gaussian noise, i.e. s{t) = n{t). Given 
the output for each configuration, we then search for the 
existence of a GW source. Even though a SMBHB source 
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does not exist in the output, the power spectrum of the 
template can match various peaks in the noise spectrum 
(and in reality, the galactic foreground of compact bi¬ 
naries also), leading to a false-positive detection. Our 
null-tests return detection thresholds of p = 9 for conhg- 
uration C2, and a threshold of p = 8 for configuration 
C3, where the two-channel SNR is defined by 


{hi\hi) 


{S2\h2f 

(^21^-2) 
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where the prior probability is given by 7r(A'^), and 
p(s|A”) = >C(A'^) is the likelihood function defined by 

£(A«) = exp (s - h{X^) \s - h{X ^))) . (12) 

Given the posterior distribution, we can now calculate 
the Bayesian credible interval (BCI), C, such that 

J p{X^\s)dX'^ = l-a, (13) 


We should mention here that, in general, the inclusion of 
a second channel usually increases the SNR by a factor 
of '/2 over a single channel interferometer of the same 
size. By comparison, as p oc L, increasing the size of 
the detector by a factor of 2 also increases the SNR by a 
factor of 2. 

To investigate the detection horizons for each config¬ 
uration, we generated populations of SMBHBs from a 
Monte Carlo simulation based on an astrophysical source 
model [28]. As with PC15, the Monte Carlo simulation 
was restricted to sources with mass ratios of g G [1,30], 
maximum array sizes of 2^^, and a time of coalescence of 
tc G [0.3,1] yr. All angular variables were examined over 
their full natural range. In Fig (2) we plot the SMBHB 
redshift detection horizon for the three configurations as 
a function of total redshifted mass. In this image, config¬ 
uration Cl is represented by the blue solid line, C2 by the 
red dashed line, and C3 by the orange dot-dashed line. 
We can see that there is only a slight improvement in the 
detection horizon for low mass systems, regardless of the 
conhguration. This is not surprising as this is the mass 
range that is dominated by the high frequency photon 
shot noise that is common for all configurations. If we 
focus on the mass range between 10^ < m{z)/MQ < 10®, 
we can see that there is an extension in the horizon max¬ 
imum from z ~ 13 for Cl, to z ~ 16 for C2 and to 
z ~ 18 for C3 (assuming that we are only using inspiral 
waveforms). These redshifts correspond to maximum de¬ 
tection limits of (141,178,203) Gpc respectively. Finally, 
we also gain in horizon distance for SMBHBs with masses 
of > 10^ Mq. This is due to both the increase in sensitiv¬ 
ity at low frequencies, and the fact that higher harmonics 
from more massive systems now contain enough SNR to 
be visible in the detector. 


IV. RESULTS 


A. Bayesian inference of SMBHBs. 


In Bayesian inference, the posterior density p(A'^|s), 
given a data set s{t), a theoretical model dependent on 
a parameter set A'^, and a noise model for Sn{f ), can be 
found via Bayes theorem 


( 11 ) 


where for a 95% BCI, a = 0.05. The BCI allows us to 
make a degree-of-belief statement that the probability of 
the true parameter value lying within the credible inter¬ 
val is 95%, i.e. 


P(A^,„, gC|s) = 0.95. (14) 

For the sky location of the source, we can define an error 
box in the sky according to [24] 


AH = (15) 

where 

= (Acos6»Acos6»), (16) 

= {A(j)A(j)), (17) 

= (Acos6»A(/)), (18) 

and = (AA'^AA'^) are elements of the variance- 

covariance matrix, found directly from the chains them¬ 
selves. 

To make a direct comparison with PC 15, we use 
the same 120 SMBHB sources out to a redshift of 
z ~ 13. In each case a combined Hessian-Differential 
Evolution Markov Chain (DEMC) was run for 10® it¬ 
erations, with a “burn-in” phase of 2 x lO"* itera¬ 
tions. During the burn-in phase, a combination of ther- 
mostated and simulated annealing was used to acceler¬ 
ate the chain mixing and convergence to the global solu¬ 
tion [29]. We also used a composite integral method [30] 
to speed up the calculations of the likelihood and the 
Fisher information matrix in the chains. For each 
DEMC, we used flat priors were used for the param¬ 
eters {InDi, InMc, In^, Infc, cos 6, COS0, (/), (/9c}, where 

the boundaries were set at Dl G [7.7 x 10“‘*,300] Gpc, 
Me G [435,5.06 X 10^] Mq, p G [250,2.9 x 10^]] Mq, 
tc G [0.2,1.1] years, {cos6, cosdj G [—1,1], (p G [0,27r] 
and were left completely open for {^p,^pc}■ The prior in 
Dl represents an assumption that there are no SMBHBs 
closer than the M31 galaxy, and no further than a redshift 
of z ^ 25. The priors in {Mc,/i} are chosen so that the 
minimum total redshifted mass is m(z) = 10® Mq, while 
the upper limit of to(z) = 1.163 x 10® M 0 assures that 
the maximum 2PN harmonic last stable orbit frequency 
corresponds to a value of 5 x 10“® Hz, corresponding to 
a prior range in reduced mass of 1 < g < 100. 


p(A'"ls) oc7r(A'=)/:(A'"), 
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FIG. 3: A comparison of 95% credible intervals for Dl (top left), Afl (top right), Me (middle left), jj, (middle right), l 
(bottom left) and te (bottom right) between a Cl (blue) and C2 (orange) eLISA-like observatory. In the cell representing the 
luminosity distance, the dashed black line represents the true injected values. For ease of interpretation, the mass parameters 
are represented as credible intervals in percentage error, while the credible intervals for (t, te) are centred around the median 
subtracted true values, and are expressed in radians and seconds respectively. 
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B. Comparison between configurations Cl and C2 

Our first objective is to answer the question of what 
do we gain by the inclusion of a third arm. In Fig. (3) 
we plot the 95% BCIs for the SMBHB parameters 
(Dl, Ar2,Mc,/x, for both the Cl (blue) and C2 (or¬ 
ange) configurations. 

Progressing parameter by parameter, we see that for 
the luminosity distance (top-left), the inclusion of the 
third arm not only reduces the size of the BCIs, but also 
extends the minimum distance that a source can be as¬ 
sumed to be at. For example, for the source at z = 4.4, 
the lower limit of the BCI for the Cl configuration is a 
distance of = 9.3 Gpc. For the C2 configuration, this 
increases to Dl = 29.4 Gpc. Similarly, at higher red- 
shifts, for the source at z = 13.2, the lower limit for the 
Cl configuration is = 23.8 Gpc, while for C2 it is 61.6 
Gpc. In general, while also reducing the upper limits of 
the BCIs, we improve the minimum possible distances of 
the source by factors of 2-3. 

For the sky error box (top-right), the third arm re¬ 
duces the size of the error box by factors of 2-4 between 
the Cl (circles) and C2 (squares) configurations. How¬ 
ever, we should highlight that depending on the position 
and orientation of the source, we can sometimes observe 
and improvement of an order of magnitude or more. As 
an example, if we take the source at z = 1.1, with the Cl 
configuration, this source has an error of AH = 594 deg^. 
The inclusion of a third arm reduces this error to 8deg^. 
We stress again that the order of magnitude improve¬ 
ments are very much source specific, and are not a general 
rule. From the scatter plot in Fig. (3), we can see that 
for this selection of sources, the binaries beyond 2 ; = 1 all 
have error boxes of > 10 deg^ for Cl, while this level of 
sky error is still attainable at 2 : = 2 for C2. At the high z 
end, while the error boxes are still large, we still acquire 
improvements of factors of ~ 2 with a third arm. 

It was already demonstrated in PC 15 that the errors in 
the other four parameters were already at a more than ac¬ 
ceptable level. However, for completion, we also include 
the results of the three arm detector here. In general, as 
expected, we see an improvement in the size of the BCIs. 
For the Cl configuration, the largest BCIs are ±1% for 
Me, ^ 18 % for /X, I 3400 secs for G and ''^^s for t. 
In contrast, for the C2 configuration, these numbers re¬ 
duce to to' 45 % for Me, for /x, ^2300 ®oos for % and 

±0.75 rads for t. 


C. Comparison between configurations Cl and C3 

We now investigate what happens if we remain with the 
eLISA configuration (blue), but double the arm lengths 
to 2Gm (orange). In Fig. (4) we again plot the 95% 
BCIs for the two configurations. The first thing we no¬ 
tice is that while the BCIs for Dl are again smaller, there 
are not as reduced as with the C2 configuration. Using 
the previous examples as a point of comparison, for the 


source at 2 : = 4.4, the minimum possible distance in¬ 
creases from the Cl value of Dl = 9.3 Gpc, to a value 
of 25.9 Gpc (as compared to the C2 value of 29.4 Gpc). 
For the source at 2 : = 13.2, the Cl observatory gave us 
a minimal possible distance to the source of 23.8 Gpc. 
This now increases to Dl = 45.6 Gpc for C3, as opposed 
to the C2 value of 61.6 Gpc. 

Now focusing on the sky error box, we can see that 
the scatter plot is more diffuse in the case of C3. This 
is to be expected as we do not have the benefit of the 
third arm. However, the reduction in the size of the 
error box is still quite significant. Again, using previous 
comparisons, the source at z = 1.1 now has an error box 
of AH = 18.6 deg^, as compared with the 594 deg^ for 
Cl and 8deg‘^ for C2. Furthermore, as with C2, the C3 
configuration also constrains the sky error to < 10^ deg^ 
out to a redshift of z ~ 2. 

As was pointed out in PC 15, due to high correlations 
between parameters, that it is impossible to talk about 
specific parameters when discussing detector configura¬ 
tions, without taking all parameters into account simul¬ 
taneously. With this in mind, while the BCIs for Dl are 
larger, and the scatter plot for AH is more diffuse, the 
BCIs for the other four parameters are in general smaller 
in the case of C3. In this case, the largest BCIs now be¬ 
come lo’ 3 ^% for Me, ±7.5% for /x, I 2100 and 

±1 rads for l. 


D. So which configuration is better? 

Looking at the results of the Bayesian inference, it is 
clear from the size of the BCIs and the sky error box, 
that both alternative configurations are better than the 
original eLISA configuration. C2 seems to give a slightly 
smaller error in both the distance and sky location than 
C3, but C3 has smaller BCIs for (Me, /x, 6, tc). Therefore, 
it seems very difficult just from an investigation of the 
BCIs alone to say which configuration is better. 

To try and shed further light on the matter, we ran a 
regression analysis based on the recovered median pa¬ 
rameter values for the 120 sources. As the errors in 
{Me, g,, i,tc) are already more than adequate, regardless 
of whether we choose C2 or C3, we will only concentrate 
on the evolution of errors in Dl and AH. In Fig. (5) 
we plot the results for the growth in percentage error for 
Dl (left) and the size of the error box (right) for the 
Cl (blue-solid), C2 (blue-dashed) and C3 (orange-solid) 
configurations. As we demonstrated above that any ad¬ 
vanced configuration would also result in an increased 
detection horizon, we plot the error growth to a redshift 
of z = 20. 

For the luminosity distance, we are usually interested 
in the redshifts below which we can usually achieve er¬ 
rors of less than 1 and 10% in the estimated parameters. 
Using this, for the Cl configuration, we achieve an error 
of 1% at z = 0.24, and an error of 10% at z = 4.23. For 
C2, these values become z = 0.51 for the 1% error, and 
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FIG. 4: A comparison of 95% credible intervals for Dl (top left), Afl (top right), Me (middle left), jj, (middle right), l 
(bottom left) and te (bottom right) between a Cl (blue) and C3 (orange) eLISA-like observatory. In the cell representing the 
luminosity distance, the dashed black line represents the true injected values. For ease of interpretation, the mass parameters 
are represented as credible intervals in percentage error, while the credible intervals for (t, te) are centred around the median 
subtracted true values, and are expressed in radians and seconds respectively. 








































z = 12.15 for an error of 10%. Finally, for C3, we acquire 
a 1% error at z = 0.16 and a 10% error also at 2 ; = 12.15. 
At z = 20, we observe percentage errors of 34.4, 14.4 and 
13.1% for the Cl, C2 and C3 configurations respectively. 
In terms of sky location error, we are again usually inter¬ 
ested in the redshift limits at which we can observe sky 
errors of AH < 1,10 and 10“^ deg^. For Cl, these red- 
shifts are z = 0.03,0.12 and 0.51 respectively. For C2, 
they become z = 0.13,0.42 and 1.33. And finally for C3, 
z = 0.11,0.38 and 1.31. We should point out that while 
C2 performs better than C3 at lower z, at z = 12 for Dl, 
and z = 2 for AH, C3 begins to outperform C2. This sug¬ 
gests that if we are more interested in distant SMBHBs, 
the C3 configuration is possibly a better choice. 

One final inference we can make is the following : if 
we assume that the inclusion of a third arm provides a 
similar parameter error improvement scaling in all cases, 
given the data from the Cl, C2 and C3 configurations, we 
can infer the performance level of a 2Gm, 6-link detector 
(which we will designate C4 in our categorization and 
is represented by the orange-dashed line in the figure). 
Using the same criteria as before, we would now acquire 
a 1% error in Dl at a redshift of z = 0.54. Beyond this, 
this growth in percentage error is so flat, that by z = 20 
the percentage error has grown to only 5.45%. For the 
sky error box, we now attain an error of AU < 1 deg^ 
at z = 0.33, < 10 deg"^ at z = 0.92 and < 10^ at 
z = 2.63. 


V. CONCLUSION 

ESA has chosen the “Gravitational Wave Universe” as 
its theme for the L3 Gosmic Vision program. Within 
this program, a mission concept called eLISA has been 
presented. This observatory will consist of one mother, 
and two daughter spacecraft, in an equilateral triangular 
configuration. The separation between each spacecraft 
will be IGm. However, while this configuration was ac¬ 
cepted for the initial proposal, the final configuration for 
the mission will not be fixed until the end of this decade. 
This allows the gravitational waves community time to 
explore alternative configurations to the eLISA mission. 

In this work, we have carried out a Bayesian inference 
analysis for 120 supermassive black hole binaries using 
two alternative configurations to the eLISA mission con¬ 
cept. In the first case we assume a three arm interferom¬ 
eter with arm lengths of IGm, while in the second case 
we assume a two arm interferometer with arm lengths of 
2Gm. In terms of scientific performance, both alternative 
configurations yield similar results and thus make it very 
difficult to assert that one is clearly better than the other. 
While the smaller three arm detector gives slightly bet¬ 


ter results for the luminosity distance and sky location, 
the larger two arm detector gives slightly smaller errors 
in the chirp mass, reduced mass, time to coalescence and 
inclination angle. 

Using the recovered median parameter values, we ran a 
regression analysis which further demonstrated the simi¬ 
larity in performance levels for the two alternative config¬ 
urations. Again, at lower redshifts, the three arm detec¬ 
tor worked slightly better than its two arm counterpart. 
But at higher redshifts this trend reversed. In both cases 
we accumulate a general error of 10% in luminosity dis¬ 
tance at a redshift of z ~ 12, growing to errors of 13—14% 
at z = 20. This is compared to the original eLISA config¬ 
uration which acquires a 10% error in distance at z ~ 4 
and an error of 34% at z = 20. In terms of sky error, we 
postulate that it should be possible to achieve errors of 
< 1 deg^ out to z ~ 0.1 and < 10 deg^ to z ~ 0.4 with 
the alternative configurations, as again opposed to the 
eLISA values of z ~ 0.03 and z ~ 0.1. 

Finally, and based on the assumption that the recov¬ 
ered median parameter values improve with a similar 
scaling factor when a third arm is introduced, we in¬ 
ferred the performance levels of a 2Gm, 3 arm detector. 
In this case, a 1% error in luminosity distance is acquired 
at z ~ 0.5, growing to an error of 5% at z = 20. We also 
improve the size of the sky error box to AU < 1 deg"^ at 
z ~ 0.3, and < lOdeg^ at z = 0.92. 

It is clear from our analysis that it is difficult to say 
whether configuration G2 or G3 is the better choice. 
However, it is promising to see that increasing the size of 
the detector is a viable alternative, in terms of scientific 
performance, to the inclusion of a third arm in a smaller 
detector. Having a third arm is always beneficial when it 
comes to localization of the source in the sky. This could 
be seen in the analysis of the Markov chains, as in gen¬ 
eral, the chains did not change modes between the true 
position and the antipodal sky position (a feature that 
was much more common in the two arm detectors). How¬ 
ever, as space based GW observatories are not telescopes, 
we expect the sky error box to increase quite quickly be¬ 
yond a redshift of z = 2. This implies that if one of our 
goals is the exploration of black hole seed formation at 
high redshift, the inclusion of a third arm may not be 
so important, especially as we have demonstrated that a 
bigger detector performs marginally better in this regime 
anyway. 

It is clear that the choice of final configurations will 
also strongly depend on technological and financial con¬ 
siderations. However, for now, the main conclusion of 
this study is that both configurations are viable alter¬ 
natives, with no significant loss of science potential for 
SMBHBs between either choice. 
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